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Abstract 

The ground-state phase diagram of two-dimensional Falicov-Kimball model with nearest- 
neighbour and next-nearest-neighbour hoppings (characterized by t, t' constants, respectively) has 
been studied in perturbative regime, i.e. in the case when on-site Coulomb interaction constant 
U is much larger than t, t' ones. The fourth-order phase diagram exhibits a rich structure, more 
complicated than for the ordinary Falicov-Kimball model. Possible experimental implications of 
the presence of nnn term are shortly discussed. 

PACS numbers: 71.10.Fd, 71.21.+a, 75.10.Hk, 75.30.Kz 
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I. INTRODUCTION 



The Falicov-Kimball model has been proposed in 1969 to description the metal-insulator 
transition lL Later on, it has been applied in another important problems: mixed valence 
phenomena crystallization and alloy formation [^] and others. In this model, we are 
dealing with two types of particles defined on a <i-dimensional simple cubic lattice Z d : im- 
mobile "ions" and itinerant spinless "electrons". There exist also other interpretations of 
the model 

The Hamiltonian defined on a finite subset A of Z d has the form 

H A = if ,A + V A , (1) 

where 

#0,A = U 2j W x Tl x - fXi 2j Wx ~ yu e ^2 Hx ' ( 2 ) 

xeA ^eA xeA 

^a = - E t{clc y + clc x ) (3) 

(x,y)eA 

Here cj. and c x are creation and annihilation operators of an electron at lattice site x G A, 
satisfying ordinary anticommutation relations. The corresponding number particle operator 
is n x = c x c x . w x is a classical variable taking values or 1; it measures the number of ions at 
lattice site x. The chemical potentials of the ions and electrons are \ii and /x e , respectively. 

The Falicov-Kimball model in its basic, "backbone" form given by (0), © is too oversim- 
plified to give quantitative predictions in real experiments. However, it is nontrivial lattice 
model of correlated electrons and captures many aspects of behaviour of such systems. It 
allows rigorous analysis in many situations; for a review, see jj] . One can hope that a good 
understanding of this simpler model might lead to better insight into the Hubbard model, 
where rigorous results are rare Jsj]. 

One can try to make the FK model more realistic by adding various terms to the "back- 
bone" hamiltonian (J2J) , Q in the manner analogous to that in the original Hubbard paper . 
(Other possibility is enlargement of the space of internal degrees of freedom 0J], but we will 
not consider it here.) The most important among them are: consideration of another types 
of lattice, particle statistics and presence of magnetic field pj; correlated hopping (analysed 
in [sj], 0, [3); taking into account the Coulomb interactions between nearest neighbours, 
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as well as (small) hopping of heavy particles ^j], j^]; consideration of the next-nearest- 
neighbour hoppings (let's name this modification as the t — if model in analogy with the 
corresponding version of the Hubbard model This last effect has been analysed in only 

few papers. In [13], authors established that if if <C t, then the phase diagram of the t — if 
FKM does not differ too much from the diagram of the pure FK model. A remarkable paper 
is devoted to analysis of three-dimensional strongly asymmetric Hubbard model (i.e. 
generalized FK one) with three hopping parameters, for large Coulomb interaction constant 
U, in the neighbourhood of the symmetry point. Authors have determined rigorously the 
structure of ground states and proved their thermal stability up to terms proportional to 
U~ 1 x (square of the hopping constant). 

In this paper, author examined influence of further terms of perturbation expansion (3-rd 
an 4-th ones) on the ground-state phase diagram in two-dimensional situation in the half- 
filling case, i.e. when the average value of the total particle number ^2 xeA (n x + w x ) is equal 
to the number of sites |A|. Effects of higher-order-terms turned out to be very interesting 
in the ordinary FK model Q], jlfi| . 

As a first step of the study, the effective Hamiltonian has been derived; it can be written 
as the Hamiltonian for the Ising model with complicated interactions, leading to strong 
frustration. After that, ground states of the Hamiltonian have been looked for, and the 
phase diagram has been constructed. In the orders 2 and 3 it was possible to determine 
it rigorously, whereas in the fourth order use of the restricted phase diagram method was 
necessary. This method proved their utility in the analysis of another versions of the FKM 
0, [3]; Q|- I n this method, one constructs collection of all periodic arrangements of 
ions up to certain values of lattice sites N per elementary cell, and then one looks for the 
configuration of minimal energy among members of this set. 

As one could expect, the fourth-order phase diagram turned out to be more complicated 
than in the case of the ordinary FKM. In this last case, one observes five phases with period 
not exceeding N = 5 j^, [3], in the case of the correlated- hopping FKM, six phases 
are present [ll3]. In the ground state phase diagram of t — t'-FKM, thirteen phases 
has been found (three of them are degenerate); the period of non-degenerate phases does 
not exceed 12 sites per elementary cell. One observes also that the region occupied by one 
phase (FK-like one, of density 1/4) is anomally large (one can expect that this region should 
occupy region of the size of the order t 4 /U 3 , whereas the actual size is of the order t' 2 /U). 
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This phenomenon can be explained by the lifting of the macroscopic degeneracy present in 
the second order by higher-order perturbation. 

Comparing phase diagrams of the FKM and t — t! FKM, it turned out that the influence 
of the nnn hopping is surprisingly large. 

The outline of the paper is as follows. In the Sec. |H1 the effective Hamiltonians up to 
fourth order perturbation theory have been derived. In the Sec. IIIIl ground states and phase 
diagram of the effective Hamiltonians in subsequent orders have been determined. Moreover, 
effects of neglected higher-order-terms as well as temperature have been discussed. The last 
section HYI contains summary and conclusions. 



II. PERTURBATION THEORY AND EFFECTIVE HAMILTON! AN 
A. Nonperturbed Hamiltonian, their ground states and phase diagram 

In this paper we examine the model in the range of parameters t, t! <C U. The value of t' 
is usually smaller than that of t, however both these quantities are of the same order. 

For derivation of the effective Hamiltonian, the method worked out in the paper 
has been applied. It has this advantage that it can serve (provided certain conditions are 
fulfilled) as a first step to application of the quantum Pirogov-Sinai method and proving 
thermal and quantum stability of ground states. Detailed description of all these procedures 
can be found in - (2^]. Here, only the application of the method and results will be 
given, as the general scheme is identical as in the paper 2^]. 

To obtain the final expression, we must divide states of the system onto ground and 
excited ones, and to find corresponding projections onto both groups. These collection of 



states are identical as in 
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Let us begin our analysis starting from the classical part of the Hamiltonian (|2*|): it 
is well known, see [22]. The Hilbert space Tii on the z-th site is spanned by the states: 
\rii t+ , Tii-) or, explicitely, |0,0), |1,0), |0, 1) and |1,1). The corresponding energies are: 
0; — — U — fi + — The phase diagram consist of the following four regions. In 
region J, defined by //+ < 0,/x_ < 0, all sites are empty. In two twin regions 
given by conditions: 27+: /i+ > 0, /i+ > //_, //_ < U (for 27_, one should interchange the 
subscripts + and — ) all sites are in the |1, 0) (corresp. |0, 1)) state. In the region 222, given 
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FIG. 1: Phase diagram of the nonperturbed Hamiltonian iffl). 

by: /i + > U, [A- > U, all sites are doubly occupied. This situation is illustrated on Fig. [TJ 

We choose the states |1,0) and |0, 1) as ground states. They are divided from excited 
ones by energy gap A = mm(U, //_) . It means that we analyse the phase diagram in 
some subset of the region II + U //_ (the shaded region on the Fig. 1). The most interesting 
situation takes place in the neighbourhood of the /i + = //_ line between regions II + and 
on this line, we observe a macroscopic degeneracy. 
The projection operator on ground states at i-th site is 



(rii 



m 



(4) 



B. Effective Hamiltonians up to 4-th order of perturbation theory 



Expression for effective 
FK model can be found in 



iamiltonian in fourth-order perturbation theory for the ordinary 
23, Table 2. The 4-th order effective Hamiltonian fort-tf FKM 
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can be derived using the same methodology, described in |22[, Sec. 3. (It should be stressed 
that expressions up to 4-th order have been derived, for the 3d model, in the paper jl^ . 
Unfortunately, authors didn't analyse effects of orders 3 and 4). 

As a final result of calculations, one obtains, after specialization to the half-filled case 
(i.e. the situation where the total number of particles is equal to the number of lattice sites): 

• Second-order correction: 

flf = -/»5> + E ^( 4 ^-!)+ E |^(4^| - 1) (5) 

i d(i,j)=l d{i,j)=V2 

where: h = fii — /i e ; Si is the classical one-half spin on the lattice site i; it is related to the 
variable Wi by the formula: Si = (wi — l)/2. 

• Third-order correction: 



^o 3) = % E ( 6 s & 3 " - 1& + s > + 



(6) 

where summation is performed over all triples {i,j,k} on the lattice such that {i,j} and 
{j, k} are nearest neighbour bonds forming the angle ir/2. 

• Fourth-order correction is the most complicated one and is a sum of two-body (2b) and 
four-body (4b) interactions: 

rr(4) _ n , f 2b \^ c-3 c-3 , T2b \^ q3 q3 , T2b \^ q3q3 , j2b V - "* c3 c3 , j2b \^ o3 c3 

d(i,j)=l d(i,j)=V2 d(i,j)=2 d(i,j)=S d(i,j)=V& 

+Ji h E/ Si^jSkSi + </ 4b E] SiSjSkSi + Jf' E] ^iSjSkSi + Jf 3 E] SiSjSkSi (7) 

ni(ijkl) ■ni(ijkV) 7T3(ijkl) iT4(ijkl) 

In formulas above, we have: C = 3r/2 + 5r' + 3r"/2; Jf h = -18r - 32r'; J 2 2b = 6r - 
36r' - 18r"; Jf b = 4r - 4r' + 6r"; J| b = 12r'; J| b = 4r"; Jf 5 = 40r + 80r'; J| b = 40r"; 
J|b = 40r ' ; J|b = 4 0r ' ; w h e re: r = t 4 /[/ 3 ; r' = t 2 t' 2 /U 3 ; r" = t' 4 /U 3 . Sets n a (ijkl) are 
defined in the following way: 7Ti ( "square" ) is formed by spins occupying vertices (0, 0), (0, 1), 
(1,1) and (1,0); ir 2 ("diagonal square") is formed by: (0,0), (1,1), (0,2) and (—1,1); 7r 3 
("big triangle"): (0,0), (0,1), (0,2) and (1,1); 7r 4 : ("rhomb") (0,0), (1,0), (2,1) and (1,1). 
The summation over four-body interactions in (J7|) is performed over all sets obtained from 
plaquettes 7Ti, . . . , 7r 4 above by operations compatible with lattice symmetries (translations, 
rotations by multiple of tt/4, reflections, inversions); the plaquette ■n a {ijkl) occupies sites 
i, j, k and I. 
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III. GROUND STATE PHASE DIAGRAMS IN ORDER 2 AND 3 



A. Order 2 

The ground-state phase diagram of the system described by the Hamiltonian (jSJ) can be 
obtained by rewriting the Hamiltonian in the following equivalent form: 



[20j |. [23[; this definition is also reminded in the Appendix). In such a case, we can replace 
the process of minimization of energy over the whole lattice by the problem much simpler: 
the minimization of energy over the set of plaquette configurations. These configurations 
are presented on Fig. 12 It leads to the picture of the phase diagram as illustrated on 
Fig. El This diagram possess two obvious symmetries. One of them is due to symmetry of 
the hamiltonian (j^J) with respect to the change of sign t' — > —t'; the phase diagram is also 
symmetric with respect to such a change of sign. The second one is the symmetry of phase 
diagram with respect to the change h — > — h; however, in this case, one should also replace 
configurations by their mirror images (i. e. S{ — > —Si). 

Phases /+ (full) and J_ (empty; for illustration, see configuration on the Fig.0J are build 
from plaquettes i + and i_ , respectively. These phases are unique. We have similar situation 
for regions III X (Neel phase; Fig. HJ configuration 10) and HI\\ ( Fig. HI configuration 11 
(it is an analogon of the "planar" phase in They are build from plaquettes iii x and 

respectively. Again, these phases are unique (modulo translations). The situation for 
phases II + , II- differs from previous ones. These phases are build from plaquettes 
ii-, respectively (see Fig. configuration No. 5 as an example); however, they are non- 
unique and possess macroscopic degeneration. One can easily check that there is a large dose 
of freedom in building of lattice configurations from these plaquettes. We encounter here 
situation similar to this which happens for the antiferromagnetic Ising model on triangular 
lattice. 




(8) 



TTl (ijkl) 



where 



= Jj( s i s 3 + SjS k + S k S t + SA) + ^(SA + SjSi) - j(Si + Sj + S k + S t ) (9) 
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I. 



II. 



Ill, 



III, 



• • • 

• • • • 

i + ii + 

FIG. 2: Possible configurations of 2 x 2 plaquettes. For plaquettes ii and Hi one should take into 
account plaquettes obtained from those illustrated above by rotations. Small dots denote empty 
lattice sites (or spins "down" in the spin language), big dots - occupied lattice sites (spins "up", 
respectively). 

B. Order 3 

Now, let us check how the phase diagram will change under switching third-order terms 
on. Let us rewrite the third-order correction (JHJ) in the equivalent "plaquette" form: 

= £ C do) 

7T1 (ijkl) 

where 

6 t 2 t' 3 t 2 t' 

^0;ijkl = jj2 (SiSjSk + SjSkSl + SkSlSi + SlSiSj) QJP^* ^ + + ^l) (11) 

(again spins k, I are arranged anticlockwise on the plaquette). 

The full Hamiltonian, up to third-order terms, is the sum of terms © and (|llj). As in 
previous Subsection, one can check that it is an m-potential. Moreover, it turns out that the 
presence of third-order terms does not modify ground states of plaquettes. In the other words, 
plaquette configurations which were ground-states in the second order, remain ground states 
also in third order! This implies that degeneracy of phases II + i iT_ is not lifted and they 
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A 

h/4 




FIG. 3: Ground state phase diagram in the second order perturbation theory for the Hamiltonian 
(jHJ), ©• h is the difference of chemical potentials. Phases I and /// are unique (modulo rotations 
and translations), whereas phases II are macroscopically degenerate. In order 3, the phase diagram 
is a small deformation of the above picture. 

still are degenerate. What does change, it is location of the boundary between phases. The 
difference in location of phase boundaries in orders 2 and 3 is of the order t 2 t' /U 2 . 

The phase diagram in third order possess certain kind of symmetry. It is discussed in 
more details in the next Subsection. At this moment, we only conclude that the phase 
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diagram in 3. order is a small deformation of the second-order phase diagram. 
C. Phase diagram in fourth order 

Regions occupied by phases II± in both second and third order exhibit macroscopic 
degeneracy. One can expect that they will be sensitive against perturbations and that in 
some of next orders this degeneracy will be lifted. This happens yet in fourth order; we 
describe the situation below. 

This picture has been obtained by the restricted phase diagram method. Recall that in 
this method, one takes into account all periodic configurations up to certain values of lattice 
sites N per elementary cell, and then one minimizes the energy over this set of configurations. 

The ground-state phase diagram is much more complicated than for the ordinary FKM - 
see Figs.|3]and|I3but it is still manageable (in some respects it is similar to the phase diagram 
of the FKM on triangular lattice, studied in [If). Thirteen phases have been detected on 
phase diagram; it seems that three of them are degenerate. Moreover, for each phase present 
for h > 0, there corresponds their "mirror" for h < (this mirror is obtained by the change 
of occupied sites onto unoccupied ones and vice versa). Periods of these phases do not 
exceed 12 sites per elementary cell. Such a picture emerged at N = 12 and hasn't changed 
up to N = 27 (which corresponds to more than 3 x 10 7 trial configurations). For this reason, 
author claims that this phase diagram is "exact but non-rigorous" . 

The phase diagram in fourth order possess certain (pseudo) symmetries. Let us describe 
the situation order-by-order. 

In the second order, we have the Hamiltonian symmetric with respect to the change 
h — > —h and t! — > —if . It is obvious that the phase diagram is symmetric with respect to 
the change h — > — h, if one change also configuration to its mirror image. 

In the third order, the Hamiltonian is no longer symmetric. However, their ground states 
are easily determined, because the Hamiltonian is expressible as a sum of m-potentials, 
defined on 2 x 2 plaquettes. It turns out that on the phase diagram, there are present the 
same phases as in the second order. The difference between phase diagrams in the second 
and third order is apparent in location of phase boundaries; corresponding lines are shifted 
by a factor proportional to tt' 2 /U 2 . This shift is symmetric with respect to the change 
h — ► —h. In the other words: If some boundary between phases i and j is shifted by e for 
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N°= 0, p = N°=l, p = 1/9 N°=l ', p = 1/9 N°=2, p = 1/8 



N°=3, p = 1/6 N°=3',p = 1/6 N°=4, p = 1/5 N°=5, p = 1/4 



• ■ • • • • ■ • 



• ■ • ■ ■ • 



N°=6, p = 2/7 N°=7, p = 1/3 N°=8, p = 1/3 N°=9, p = 2/5 



N°=10, p = 1/2 N°=U, p = 1/2 N°=12, p = 3/8 N°=12', p = 3/8 

FIG. 4: Phases appearing on the ground-state phase diagram for the fourth-order effective Hamilto- 
nian Q. Configurations: 0,2,4-11 are unique modulo lattice symmetry operations. Configurations 
possessing densities 1/9, 1/6 and 3/8 are degenerate and form (infinite?) series; two first examples 
for every such a series are shown (1,1' etc.). 

h > 0, then, for h < 0, the boundary between mirrors of phases i and j is shifted by — e. 

Let us analyse (pseudo) symmetries of the fourth order phase diagram. Consider first 
some third-order configuration of minimal energy for h > (call it CME^); for every such a 
configuration, we have its mirror, which is also the configuration of minimal energy for some 
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A 

t'/t 




8t 2 /U h 



FIG. 5: The ground-state phase diagram for the fourth-order effective Hamiltonian (JJJ). Only 
topological aspects of the phase diagram are displayed, as most of phases occupy very narrow 
regions. Phases: 1-4, 6-9 and 12: width of these regions is of fourth order in expansion parameter; 
other phases occupy regions of width of the second order in expansion parameter. Only one 
quarter (i.e. values h > 0, t' > 0) are shown, as for other quarters, topological structure of the 
phase diagram is the same. 

h < (call it CME?.). Now, let us add the fourth order contribution to the Hamiltonian. 
This term is symmetric with respect to reversing of spins; so, if for h > we have some 
configuration of minimal energy CME+, then for h < the configuration of minimal energy 
CMEi will be the mirror of CME^. 

Let us summarize the above by the statement that the topological structure of phase 
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diagram is the same for h > and h < 0; for every phase i appearing for h > 0, we have 
corresponding mirror i for h < 0. Phase boundaries between phases i and j and their mirrors 
i and j are related by: 

K/j = At 2 + A't' 2 + BtH' + Ct A + C'th' 2 + C'Y 4 

^ = -At 2 - A't' 2 + BtH' - Ct 4 - C"t 2 t' 2 - C"t' 4 

Let us stress that this situation (i.e. the same structure of phase diagram for h > and h > 
is a very peculiar property of the fourth-order Hamiltonian ©, ©, (J7J); it is due to degen- 
eracy and (pseudo) symmetries of lower-order Hamiltonians. For more general Hamiltonian, 
we have no such similarity. 

Perhaps, the most interesting effect of the presence of the term with next-nearest- 
neighbour hopping is the appearance of the anomally large region occupied by the phase 
5 (FK-like phase with density 1/4). At first sight, phases appearing in the fourth-order 
perturbation theory should occupy region of the width p(t,t')/U 3 (where p(t,t') is some 
homogeneous fourth-order polynom in t and t'). However, it turns out that phase 5 occupies 
a region of width proportional to t' 2 /U, i.e. of the same order as the Neel phase, appear- 
ing in the second order! One can explain this phenomenon in the following way: Regions 
occupied by phase II + (the situation with phase II- is analogous) has width of the order 
t' 2 /U proportional to exhibit macroscopic degeneracy both in 2-nd and in 3-rd order. The 
fourth-order perturbation lifts this degeneracy and as a result the phase 11+ (an 'ancestor') 
transforms into non-degenerate phase 5 (the 'descendant') of the same density, and occupies 
region approximately as large as an 'ancestor' II + . 

Most of phases is unique, but there are also phases which remain degenerate even in 
the 4-th order (phases 1, 3, 12). Strictly speaking, the restricted phase diagram method 
detects here only finite degeneracy, i.e. finitely many ground states with identical energy 
and density but different orderings (we don't count trivial degeneracy due to symmetry 
operations, i.e. translations, rotations and reflections). Number of these ground states 
grows with N; for N = 27 we observed: eight phases of equal energy and density 1/9, (the 
first two such configurations are phases 1 and 1'); eight phases of density 1/6 (the first two 
such configurations are phases 3 and 3'); five phases of density 3/8 (the first two of them are 
12 and 12'). Moreover, it has been observed that every member of such collection of phases 
with equal density is build up from identical plaquette configurations 3x3. They can be 
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'glued' together in various arrangements and there is no uniqueness in such a procedure, i.e. 
resulting lattice configuration is non-unique. Situation here is similar to that which happens 
in the second order for phases II±. It is natural to conjecture that we encounter here the 
macroscopic degeneracy, i.e. presence of infinite number of configurations of identical density 
and energy. However author (so far) can't prove this. 

Phase diagrams in orders 2 and 3 are rigorous (by writing out Hamiltonians as sums of 
m-potentials). Author tried to do analogous thin g in the fourth-order by an attempt to 
construct m-potentials in a manner analogous as in jl^. Q, 1(3] > however, it succeeded only 
for some phases, but not for a whole phase diagram. 

An analysis above concerned the 'truncated' phase diagram, i.e. the phase diagram of 
the fourth-order effective Hamiltonian (J7j). Which changes can result as an effect coming 
from next (neglected) orders and temperature? Author expects (and conjectures) that most 
phases (i.e. all with exception of degenerate ones: i.e. 1, 3, 11 and their descendants) 
are stable ones Ql 13 • For such phases, Peierls conditions (both classical and quantum 
ones - see jioj]. are fulfilled, and regions occupied by these phases deform in only 

small manner upon thermal and quantum perturbations. This assertion concerns regions of 
phase diagram sufficiently far from phase boundaries. For regions of width of the order of 
t 5 /U A around phase boundaries, one cannot formulate any statements without going into 
next orders of perturbation theory. However, at present author cannot prove stability and 
Peierls conditions. (It could be proved by construction of m-potentials, which failed so far). 
Regions occupied by degenerate phases constitute different problem. Here, changes caused 
by thermal and quantum perturbations can be very significant; they are "terra incognita" 
and we skip this subject, leaving it as an open problem. 



IV. SUMMARY AND CONCLUSIONS 



The effective Hamiltonian and phase diagram for ground states of the t — t! FKM have 
been determined up to fourth order of perturbation theory. In the second and third order, 
phase diagram was constructed by rewritting the Hamiltonian as a sum of m-potentials. The 
phase diagram in the fourth order has been determined by the method of restricted phase 
diagrams; author claims that such a picture is "exact but non-rigorous" . The phase diagram 
is considerably more complicated than for the ordinary FKM, but still it is manageable. 
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Thirteen phases are present (plus their "mirrors"); three of these phases are degenerate 
(i.e. possess identical density and energy but different ordering). Author conjectures that 
this degeneracy is macroscopic (the method of restricted phase diagrams handles only finite 
number of configurations). 

Let us list some of the features of the phase diagram in fourth order: 

1. For small t' , the phase diagram is similar to this for the ordinary FKM; however, for 
larger t', these phase diagrams are quite different. 

2. One observes anomally large region occupied by one of the phases appearing in 4-th 
order (phase number 5). It could have experimental implications as a possibility of 
appearance of "charge density waves" more exotic than Neel ordering. 

3. The phase diagram of the full model (i.e. with inclusion of neglected terms perturba- 
tion theory, as well as temperature) is more difficult to examine than for the ordinary 
FKM. At this stage of investigation, almost nothing rigorous can be said. However, 
author expects that after switching the quantum and thermal perturbations on, the 
phase diagram will change in only small manner inside regions occupied by nonde- 
generate phases (0,2,4-10,12). But regions of width of the order t 4 /U 3 around the 
phase boundaries, as well as regions occupied by degenerate phases 1,3,11, are out of 
possibilities of present analysis. 

4. The Falicov-Kimball model is sensitive to perturbations. Modifications of the Hamil- 
tonian such as introduction of correlated hopping, nnn-hopping or consideration of the 
FKM in non-perturbative regime (i.e. for values oit/U not very small) can significantly 
or even drastically modify the phase diagram. 

APPENDIX A: GROUND STATES, M-POTENTIALS 

The Hamiltonian , where A C Z d , is a function defined on Q,\x\ - the space of all 
configurations of the system. Usually Hamiltonian is defined as a sum of potentials, i.e. 
functions defined on subsets of A: = J2bca^ b - Usually one imposes restrictions such 
that potentials are finite-range ones, i.e. such that \B\ < M, M finite. It is also assumed 
that potentials are translation invariant. 
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m-potential. Now, consider the system on an infinite lattice Z d . Assume that sets 
B a (i.e. potential supports) are translation of a fixed plaquette B by a lattice vector a: 
B a = r a B, where r a is an operator of such translation. We say that the function is an 
m-potential, if there exist configuration (perhaps, non-unique) uj° G fl^ with the following 
properties: i) For every plaquette B a , the "plaquette energies" , i.e. values of the Hamiltonian 
calculated on the plaquette B a : $B a (c<j ) are all equal; ii) For every another configuration 
uj G fijv, the condition: & Ba (u°) < &B a (u) is fulfilled. 

Ground states of the classical Hamiltonian. If there exist such configuration u° G 
Qn as above, then we call it the ground state of the Hamiltonian. 

The property of the potential to be an m-potential can be reformulated as follows. If 
a given potential is an m-potential, then the local minimality of energy (i.e. minimality 
on a plaquette) implies the global minimality on the whole lattice. The property that a 
potential is an m-potential is very important one, as it replaces searching of ground states 
of the infinite lattice by looking for the minima on a finite set. Unfortunately, some given 
potential possess this property only exceptionally. (But fortunately, in Sees. IIII Al and IIII Bl 
they share such property!) A method to avoid this obstacle is to find - for a given potential 
$ - an equivalent potential <&', such that $' is an m-potential. However, in general it is 
difficult task. 
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